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ABSTRACT 

Planets with sizes between those of Earth and Neptune divide into two populations: 
purely rocky bodies whose atmospheres contribute negligibly to their sizes, and larger 
gas-enveloped planets possessing voluminous and optically thick atmospheres. We 
show that whether a planet forms rocky or gas-enveloped depends on the solid sur¬ 
face density of its parent disk. Assembly times for rocky cores are sensitive to disk 
solid surface density. Lower surface densities spawn smaller planetary embryos; to as¬ 
semble a core of given mass, smaller embryos require more mergers between bodies 
farther apart and therefore exponentially longer formation times. Gas accretion sim¬ 
ulations yield a rule of thumb that a rocky core must be at least 2before it can 
acquire a volumetrically significant atmosphere from its parent nebula. In disks of 
low solid surface density, cores of such mass appear only after the gas disk has dissi¬ 
pated, and so remain purely rocky. Higher surface density disks breed massive cores 
more quickly, within the gas disk lifetime, and so produce gas-enveloped planets. We 
test model predictions against observations, using planet radius as an observational 
proxy for gas-to-rock content and host star metallicity as a proxy for disk solid surface 
density. Theory can explain the observation that metal-rich stars host predominantly 
gas-enveloped planets. 

Key words: planets and satellites: formation planets and satellites: terrestrial 
planets - planets and satellites: dynamical evolution and stability 


1 INTRODUCTION 

The Kepler Mission has discovered an abundance o f 
planets — of order one per star dFressin et al.l l2013l l 
- with orbital periods shorter than a year, and siz es 
smaller than that of Neptune (fBorucki et al.l 1201 1 alibi: 


Batalha et ah | 120131 : iBurke et al.l 
20151 '). Kepler measurements of 
combined with planet 


2014 iMullallv et al l 
planet radii, when 


_ asses measured via radial- 

velocity follow-u p (iMarcy et al.| l20l4 : IWeiss fe Marcvll20l4 : 

[Dressi^^^^^^^^^^3^^^3}_and_transittiming_varia^ 

tions ^ord_et_alJ_[ 20 UEj£^kr^lS r ^^iIli 221 ^i--lST^-^ . a V^ 


2012a| tL|ljbhwiclLxie^.Wu|_ 20121 SteRnmt_alJ_l2012b|la 


Mazeh_et_al .ll2013l:'l Stcffcn et al . 2013i : IWu fc Lithwickll20ia 


Had den & Li thwick 2Q14; Xiel 20141 1 , reveal a surprising di¬ 
versity of bulk densities. Some planets are dense enough 
to be made practically e xclusively of so lid refractory ele- 
ment s — e.g., Kepler-10 b llBatalha et al.|[201ll l and Kepler- 
36b (ICarter et alj 12012 1 — whereas others require volu- 
minous volatile componen ts, e.g., the Kepl er-11 planets 
dLissauer et al. 1 1 20 111 . l2013l l and Kepler-36c dCarter et aid 
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l2012ll . We refer to the former class as “purely rocky” (a.k.a. 
“super-Earths”) and the latter as “gas-enveloped” (a.k.a. 
“mini-Neptunes”). 

Purely rocky planets can undergo transformation 
to gas-enveloped, and vice versa, throug h a variety 
of physical processes incl uding outgassing f R,ogere_ek_alJ 
201l|f , photoevap oration dLopez. Fortney &: Miller] 20121 : 
Lopez fc For tnevl 20131), and collisional strip ping of atmo¬ 


spheres ( Schlichting, Sari fc Yalinewicbll2015T ). Our interest 
here is in determining a planet’s make-up just after its for¬ 
mation. We seek to identify the factors that dictate whether 
or not a rocky core can accrete a vol umetrically significant 
atmo sphere from its parent nebula dLee, Chiang fc Ormell 
l2014lh We highlight the decisive role played by the surface 
density of solids in the primordial disk — and by exten¬ 
sion the host star metallicity — in controlling how long a 
rocky core takes to coagulate. According to our interpre¬ 
tation, purely rocky planets are born within disks of low 
solid surface density and accordingly have long gestation 
times: they attain core masses large enough to acquire at¬ 
mospheres only after their ambient gas disks have dissipated. 
Recent breakthroughs l i nking planet radius to composition 
dLopez fc Fortnevl I20l4 iRogerd 1 201 5l ) allow us to identify 


© 2015 RAS 












































































2 Dawson, Chiang, & Lee 


with confidence purely rocky vs. gas-enveloped planets in 
the Kepler sample for comparison with theoretical models. 

We present a planet formation scenario in Section 2, 
deriving approximate analytic estimates for the masses of 
the seed bodies or “embryos” that merge to form super- 
Earths and the underlying cores of mini-Neptunes, and the 
ambient gas densities during the merger era. These order-of- 
magnitude considerations set the stage for numerical simu¬ 
lations in Section 3, where we combine iV-body integrations 
with atmospheric accretion models to show that disks with 
the highest solid surface densities spawn gas-enveloped plan¬ 
ets, whereas disks having sufficiently low solid surface den¬ 
sities produce purely rocky planets. In Section 4, we trans¬ 
late, as best we can, model predictions for gas mass frac¬ 
tion vs. disk surface density into the observables of planet 
radius vs. host star metallicity. There we compare to the ob¬ 
served sample of Kepler candidates with spectroscopic host 
star metallicities dBuchhave et alJ 120141 1 and find support 
for the link between disk surface density and small planets’ 
composition. We conclude in Section 5. 


2 OVERVIEW OF FORMATION SCENARIO 
AND INITIAL CONDITIONS 

The solids in a protoplanetary disk provide the raw material 
for building rocky planets and the cores of gas-enveloped 
planets. Planetesimals may have either formed in situ or 
accumulated by radial drift from larger orbital distances. 

In a disk of solid surface density E z , an “embryo” 
(a.k.a. isolation mass) consumes all the solids within its feed¬ 
ing zone of annular width A a: 

M em b = 27raAaE z . (1) 


To give a sense of scale, if Aa = 5Rn, where 


f 2M em b 
V 3 M* 


1/3 


( 2 ) 


is the mutual Hill radius and A'K is the host star’s mass, 
then Me m b is approximately the mass of Mars at semi-major 
axis a = 0.3 AU in a disk with E z = 30g/cm 2 (o/AU) _3//2 
(for context, this surface density is about 3x that of the 
solid component of the m inimum-mass solar nebula; see, e.g., 
IChiang fe Youdinll2010l and use their Z le i = 0.33). We will 
consider other possibilities for A a and E z below. 

At short orbital periods, isolation-mass embryos grow 
quickly, consuming the planetesimals in their annuli well 
before the gas disk dissipates. The growth or coagulation 
timescale is approximately 

, Pe mb Remb ( , ' / 0 , 

tcoag ~ N 

where Q is the orbital frequency, p em b the embryo’s den¬ 
sity, f?emb the embryo’s radius, and JF > 1 is the enhance- 
ment of the growth rate due to g ravitational focusing (e.g., 
iGoldreich, Lithwick fc Saril 120041 ). For example, even with¬ 
out gravitational focusing [T = 1), the Mars-mass embryo 
mentioned above takes only 0.2 Myr to grow. 


2.1 Embryo self-stirring as the gas disk dissipates 

While the gas disk is still present, it damps the embryos’ 
eccentricities by dynamical friction, preventing them from 
crossing orbits, merging, and growing to larger masses. The 
gas disk needs to dissipate so that the embryos’ self-stirring 
becomes more effective against dynamical friction. The goal 
of this subsection is to gauge the degree of gas depletion, 
and the magnitudes of the embryos’ random velocities and 
spacings, at the time that embryos cross orbits - i.e., at 

the onset of the “giant impact” stage of planet formation. 
This will inform the choice of initial conditions for the N- 
body coagulation simulations in Section 3 — which will not 
include gas dynamical friction. 0 

We sketch below the various regimes for the embryos’ 
equilibrium random velocities v. There are three velocity 
scales. The Hill velocity, «h, tends to be the smallest: 


«H 


( 2 M emb \ 1/3 
V 3 M* ) 


0.4km/s 


(M emb \ 1/3 

V M @ ) 




( 4 ) 


The sound speed, Cs, is not much larger. We assume a disk 
temperature T of 


T = 1500K 


0.1 AU 


1/2 


max(a, 0.1AU) 


( 5 ) 


yielding a sound speed 


c s = y/kT / p = 1.29 km/s 


1 AU 


1/4 


max(a, 0.1AU) 


( 6 ) 


where k is the Boltzmann constant and 
p = 2.34 x 1.67 x 10~ 24 g is the mean molecular mass. 

The largest scale is the escape velocity from the surface 
of the embryo: 


Vesc 


2GM em h 

^emb 


11 km/s 


/ Memb \ 

\ M® J 


1/2 



( 7 ) 


1 Our analysis in Section 2 assumes that embryos do not open 
gaps in the g as disk . Whe ther this is true depends on the disk 
viscosity (e.g., Duffell2015l. and references therein). Gaps weaken 
eccentricity damping by gas, but open the po ssibility of eccentric¬ 
ity excitation by gas iGoldreich fe Sarill2003h . In any case, we do 
not expect gaps to seriously interfere with the ability of rocky 
cores to accrete gas envelopes, since the gas accretion rate is in - 
sensitive to the ambient gas density llLee, Chiang k, Ormelll2Q14b 
and since gas in the gap is continuous ly replenished by viscous 
radial inflow llLubow fc D’Angeldl200(ll . 
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Rocky Recipe 3 


2.1.1 Growth of random velocities to uh 


2.1.4 Gas damping shuts off 


During the first stage of stirring, the embryos’ random ve¬ 
locities are below uh- Their stirring rate is 

3 


) Cl,v < (8) 

uh 


stir 1 Pemb ^emb ^ 

llGoldreich. Lithwick fc Saril l2004j ). At the same time, gas 
damps random velocities at a rate 

4 

Cl, v < c s (9) 


^esc 

C s 


damp, 1 Pemb-Remb 

fe.g.. iKominami fc Ida 2002), resulting in an equilibrium, 

s, (c s y 

'iJequi,! I \ Cs^V UH • 

S ga s \ 'CH J 

The embryos’ random velocities grow to vh once the gas 
surface density drops to 

\ 4 

S ffas < s 2 [ — 

uh 


This condition evaluates to E ga s < 700E Z for Mars-mass em¬ 
bryos at a — 0.3 AU, or E ga s < 100E Z for Earth-mass em¬ 
bryos at 1 AU. For disks initialized with cosmic gas-to-solid 
ratios of E gas /E, ~ 200, little or no depletion is required for 
v ~ vh. 


2.1.2 Growth of random velocities to c s 


Next the big bodies stir each other from random velocities of 
wh to c s . The gas damping rate remains the same (Equation 
EDi. but the stirring rate decreases to 


stir,2 


Pemb Remb 




VH < V < V e: 


( 10 ) 


(|Goldreich. Lithwick fc Sarilhooj ). resulting in a new equi¬ 
librium 


/ E~ \ 1,/4 

^equi,2 ^ “ 1 C s , VH “C V <C C s . 

\ ^-/gas / 


Random velocities reach v ~ c s once E gas < E z . 


2.1.3 Growth of random velocities to %/3c s 

Once the gas density drops enough for v to exceed Cs, the gas 
damping rate decreases to the classi cal formula for dynam¬ 
ical friction ( IPapaloizou fc Larwoodll2000l ). We assume the 
random energy is equipartitioned so that the random ver¬ 
tical velocity is c s /\/3 and the embryo remains marginally 
embedded in gas. The new damping rate is 


Finally, once S gas < E z /\/3 so that v > y/Scs, the embryo 
is no longer embedded in the disk but plunges through the 
disk twice per orbit, reducing the gas damping further to 



Egas 

( Uesc ^ 

W damp,4 

Pemb -Re mb 

V V J 


G, V3 c s < 


v 


( 12 ) 


llFord fc Chiand [20071 : iReinl l2012l 'l . At this stage, stirring 
(Equation 1101) always exceeds gas dynamical friction; gas 
damping effectively shuts off and cannot establish a velocity 
equilibrium. 


2.1.5 The case of widely spaced embryos 


The stirring rates cited above (Equations 8 and 10) are based 
on the particle-in-a-box approxi mation, which holds o nly for 
embryos spaced closely enough (iFord fc Chianell2007l l. Spac- 
ings widen as embryos merge, eventually invalidating the ap¬ 
proximation. We estim ate the reduced stirring rate by using 
the encounter map of lHaseeawa fc Nakazawal dl99(j ). Two 
planets separated in semimajor axis by uRh undergo close 
encounters every synodic period A t s = 4na/(3nRnTl). Each 
encounter increases planetary eccentricities by Ae as given 
by 

A(ae/f? H ) 2 « 2(a/R H ) 2 eAe « (13) 

n 4 

(jffasegawa fc Nakazawalll99f)l . their equation 38). The stir¬ 
ring rate is then 


Substituting 


v 1 Ae 135R|f2 
v e t B 8nn 3 e 2 a 3 


E z = 


^embPemb 

3anRn 


yields 


Pemb-Remb 


■Cl. 


(14) 


(15) 


(16) 


This formula is probably an upper limit to the true stir¬ 
ring rate, as it assumes that the planets’ eccentricity and 
inclination vectors are randomized between encounters (i.e., 
Equation m is “phase-averaged”). 

For v < c s , equating Equation 1161 to the damping rate 
in Equation 9 results in a velocity equilibrium 


Uequi,5 


E gas nv H 


(17) 


f V \ E gas { Vesc \ / Vesc \ 3 q 

\ V / damp,3 Pemb^emb V C s / ' v ' 

Cs < V < V3Cs, 


( 11 ) 


and the random velocity reaches a new equilibrium of 


Uequi,3 


^-Cs,Cs < V < %/3Cs 

^gas 


Until Sgas rsJ S 2 : /V3. 


The equilibrium velocity reaches Cs when the gas surface 
density drops to 



This condition evaluates to E gas < E z for Mars-mass em¬ 
bryos at 0.3 AU separated by n = 5. 

Once the equilibrium velocity exceeds c s , the ratio of 
the stirring rate (Equation 16) to the damping rate (Equa¬ 
tion 11) grows with v, and the velocity equilibrium is no 
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4 Dawson, Chiang, & Lee 


longer stable. In this case, gas damping shuts off. However, it 
may not stay shut off. Equation 1181 suggests that as embryos 
grow more widely spaced through successive mergers, eccen¬ 
tricity damping by gas may return to significance, bringing 
the equilibrium velocity back below c a and necessitating fur¬ 
ther reductions in E gaa to effect more mergers. For example, 
Equation[l 8 ]evaluates to E gaa < 0.04E 2 for 2M® cores sepa¬ 
rated by 10 Ra at 0.3 AU. Intermittent gas damping would 
both increase the total core growth timescale and limit the 
amount of gas a planet can ultimately accrete. We will not 
account for such ongoing gas damping in the remainder of 
the paper, but a detailed, self-consistent treatment for pa¬ 
rameters appropriate to close-in sup er-Earths is an impor- 
tant subject for future work (see, e.g., iKominami fc Id~all2002i 
for a study pertinent to solar system terrestrial planets). 


2.2 Growth to super-Earths and mini-Neptunes 
in the giant impact stage 

The results of the previous section suggest that, as a rough 
rule of thumb, E gaa ~ E z is the condition that deactivates 
gas dynamical friction and that triggers when embryos can 
cross orbits and merge; furthermore, at the start of this gi¬ 
ant impact stage, random velocities v ~ c a . An a lternative 
scalin g for the random velocity is to use uh; llda &; Makind 
d 199311 find that Hill scalings are appropriate for mutual stir¬ 
ring of oligarchs in the presence of planetesimals (neglecting 
gas). We will experiment with both scalings, respectively, in 
sections 12 . 2.11 and 12 . 2.21 

With E gaa ~ E 2 , the giant impact stage commences 
with ample gas for merging embryos to accrete. How much 
gas the growing rocky cores actually accrete depends on how 
long the merging process takes. If a super-Earth-mass core 
takes too long to coagulate, the already depleted gaseous 
nebula will have completely dissipated (by viscous accretion 
onto the star) before the core can acquire a volumetrically 
significant atmosphere. 

T he timescale T for em bryos to cross orbits is 
(iYoshinaga. Kokubo fc Makind 1 19991 1: 

log 10 T = Ci + C 2 n (19) 

where C i and C 2 are constants that depend on the initial 
random velocities — in units of «h — and n is the number 
of mutual Hill radii Rh between embryos. Therefore for a 
given set of embryo separations and initial velocities scaled 
by their Hill values, the timescale for the first merger is in¬ 
dependent of E 2 . But for the small values of M em b in low E 2 
disks, one merger is not enough to build super-Earth cores; 
the timescale to grow to M core > M em b depends strongly on 
S 2 , as we now show. 

2.2.1 Initial conditions set by the Hill scale 

We start by considering the case where the initial 
Ci, C 2 , and n are in depe ndent of S 2 . For exam ple, 
lHansen fc Murravl ll2012l l and lHansen fe Murray! d2013h re¬ 
cently simulated the giant impact stage with initial orbital 
spacings defined using a fixed multiple of Rh, and random 
velocities independent of E 2 . In this set up, the first orbit 
crossing time T is independent of S 2 . However, the number 


of mergers will depend on E 2 . For an initial spacing of no 
Hill radii, the embryo mass (Equation [lj is 



where E 2 = E Zi i(a/AU) _ ° (e.g.. iKokubo fc fdall2002h . 

To grow to M core — here loosely defined as the mini¬ 
mum mass necessary to accrete a volumetrically significant 
atmosphere — equal-mass embryos must undergo approxi¬ 
mately log 2 (M core /Memb) mergers, each of which increases 
the number of Hill spacings n by 2 2 ^ 3 (since each merger 
increases the absolute separation between embryos by a fac¬ 
tor of 2 while the mutual Hill radius increases by a factor 
of 2 1 / 3 ). The timescale for the final merger to form M core is 
therefore 

log 10 2fin a l = Cl + C 2 (Mcore/M em b) 2//3 no, 
which we can rewrite as 

log 10 Tfina.1 = Ci + G 3 /E 2 , 

where C 3 is independent of E 2 . Therefore the time to form 
M core depends exponentially on the solid surface density. 


2.2.2 Initial conditions set by gas sound speed 


If the initial spacings and random velocities are scaled in¬ 
stead by the gas sound speed c a (sections 12.1.2112. 1.411 . then 


Af e mb ^ 27raE 2 C a /fl 


and 


no = 


c s 

HRh 


The timescale for the final merger is 


logiQ 


C 1 +C 2 


-Mcore Mq^ 3 

M e mb 


C s 

f la 


The coefficient C 2 may also depen d 011 S 2 , because c a doe s 
not scale with «h- IYoshinaga. Kokubo fc Makind (il999f) 
found that C 2 « 0.84(1 — 0.8(e 2 ) 1 1 / i ^ ial ), where e = ea/Rn. 
Then for e ~ c B /vn, 


log 10 1 f]nal — Cl -(- 


c 3 



where C 3 and C 4 are independent of E 2 . Alternatively, if 
initial random velocities are damped to zero by mergers that 
precede the final doubling, then 


log 10 ^flnai — Ci + C 3 /E z . 


Either way, we see that the final merger to reach M core takes 
a time that is exponentially sensitive to E z . 


3 SIMULATIONS OF CORE GROWTH AND 
GAS ACCRETION 

We saw in Section 2 that embryos begin crossing orbits 
and merging when the gas surface density E gaa drops below 
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Rocky Recipe 5 


the solid surface density E z . The subsequent merger history 
unfolds over a timescale that is exponentially sensitive to 
E z . High E z disks promise to spawn massive cores quickly 
enough that they can accrete volumetrically significant at¬ 
mospheres before the gas disk dissipates. Here we determine 
using IV-body coagulation simulations whether this promise 
can be fulfilled. 


3.1 TV-body simulations of core growth 


We run fou r sets of IV-body integrations using mercury6 
(Ichambersl 1 19991 ) to si mulate the era of gia n t imp acts, 
similar to those run bv IChambers fe WetherilJ (119981 ') and 
lHansen fc Murravl (|2012| . 2013T ). We use the hybrid symplec- 
tic integrator with a time step of 0.5 days and a close en¬ 
counter distance (which triggers a transition from the sym- 
plectic integrator to the Burlisch-Stoer integrator) of 1 Ru¬ 
in each set of 500 integrations, we begin with isolation-mass 
embryos resulting from disk solid surface densities spanning 
3-400 g/cm 2 at 1 AU (= E z p). All surface density profiles 
scale as a -3 / 2 . 

In the first set of simulations (corresponding to the 
scaling arguments in Section 2.2.1), the embryo masses 
are defined by annu li each spanning 1 0 mutual Hill radii 
(10-Rh), simila r to Hansen fc Murravl (120121 . 120131 ) and 
iKokubo fc Idal (1 199£ . 20021 1. The embryos’ initial semi¬ 
major axes span 0.04 < a < 1 AU. Initial random veloci¬ 
ties are set to the Hill velocity vh, equipartitioned between 
eccentricity (e = y/2/3vu) and inclination (i = uh/%/3). 
In the second set of simulations (corresponding to Section 
2.2.2), the embryo masses are defined by annuli spanning 
2c s /n, and random velocities are set to c s (Equation [6]), 
similarly equipartitioned. In these first and second sets of 
simulations, all bodies have bulk densities set to 1 g/cm 3 , 
which defines their collisional cross-sections. The third and 
fourth set of simulations correspond to the first and second 
sets, respectively, except with bul k densities of 5 g/cm 3 in - 
stead of 1 g/cm 3 . (For comparison, lHansen fc Murravl (|2013l ') 
use 3 g/cm 3 .) Because the results did not change materially 
upon using the higher bulk density, all of the results shown 
below will be drawn from our first and second simulation 
sets. 


The solid surface density, by means of setting the ini¬ 
tial embryo masses, strongly affects the core coagulation 
timescale. Smaller initial embryos necessitate more merg¬ 
ers to reach a given core mass. More mergers take longer 
time (exponentially longer, according to the rough argu¬ 
ments presented in Section 2). In Figure[l] we show example 
snapshots from three simulations (all with Hill-scaled initial 
conditions) of disks with different E z . Cores quickly reach 
substantial masses (M core > 2Af®) in the highest E z disk; 
take millions of years to form in a lower E z disk; and remain 
low mass in the lowest E z disk. Next we summarize the re¬ 
sults of two sets of 500 simulations a piece. In Figures [2]and 
[3] we plot the time to grow to M co re, as a function of E z ,i, 
for four sample values of M core . Figure [2] corresponds to Hill- 
scaled initial conditions and Figure [3] to those set by c B . A 
horizontal dashed line marks an approximate remaining life¬ 
time of 1 Myr for the gas disk once E gas has declined to E z . 
We employ only those simulations for which the maximum 



Figure 1 . Snapshots of coagulating protocore masses vs. semi¬ 
major axis from three simulations, color-coded by solid surface 
density normalization E z i. For reference, the minimum-mass so¬ 
lar nebula has E z ] ~ 11 g/cm 2 . Snapshot times are listed verti¬ 
cally on the right. Planets in the high solid surface density sim¬ 
ulation (orange diamonds) grow quickly and reach 2(dotted 
horizontal line) well within the gas disk lifetime of ~1 Myr, ma¬ 
triculating to become gas-enveloped mini-Neptunes (large sym¬ 
bols). For the intermediate solid surface density (blue circles), 
planets grow to 2M0 only after millions of years. In the lowest 
solid surface density disk (black squares), planets fail to reach 
2M® over the 27 Myr duration of the simulation. 
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6 Dawson, Chiang, & Lee 


Table 1 . Parameters for models of gas accretion in Figure [4] 


core mass attained is less than 30 Earth masses. For a wide 
range of E 2 , the embryos can grow to M core in sufficiently 
short time to accrete atmospheres and become enveloped 
by gas; in some cases, the embryo mass itself exceeds M core 
(gray regions). But for a subset of lower E 2 disks, the time 
to grow to M core is much longer than the remaining gas disk 
lifetime; such planets run out of time to accrete gas and 
remain purely rocky. 


3.2 Atmospheric accretion models 

iLee, Chiang fc Ormell (l2014l l computed how a rocky core of 
a given mass accretes gas from the primordial nebula, with 
specific application to Kepler planets. We use their model to 
generate, for an array of core masses, a set of gas accretion 
histories that we will combine in Section 3.3 with the N- 
body coagulation simulations from Section 3.1. This combi¬ 
nation will enable us to determine, for any given protoplanet, 
the simultaneous time evolution of core mass M cole (t) and 
atmospheric mass fraction GCR(t) = gas-to-core mass ratio. 

Fo r each core mass M co re £ {1, 2, 4,10}A7®, we run the 
iLee, Chiang fc Ormell d2014ll model for two nebular temper¬ 
atures T £ (900, 600} K appropriate for two orbital dis¬ 
tances a £ {0.3, 0.7}AU, respectively. For every T and M core , 
we adopt two nebular gas densities. The higher gas density 
p gaB = E gas n/c s is chosen such that E gas = E 2 , where E 2 is 
that (extreme) value required to form M core as an isolation- 
mass embryo from a lO-ife-wide annulus at the given a. 
The lower gas density is 1/100 this value. The set of models 
so constructed span a large enough range of p gaa and M COTe 
that they can be usefully interpolated. We run one set of 
4[M core ] x 2[T(a)] x 2[p gas ] = 16 gas accretion simulations 
using solar metallicity opacities with dust, and another set 
of 16 simulations using solar metallicity o pacities without 
dust (for technical detail s, see ILee. Chiang fc Ormell 120141 
and I Ferguson et al.l l200ol ) . Model parameters are summa¬ 
rized in Table [T] including the corresponding values of E 2 ,i. 

Figure [4] plots the resulting gas accretion histo- 
ries for dusty-atmosphere opacities. As discussed by 
ILee. Chiang fc Ormell ( 2014 1. these histories are not very 
sensitive to outer boundary conditions like T, a , or p gas , be¬ 
cause accretion is controlled by the ability of the atmosphere 
to cool internally, i.e., by conditions at the atmosphere’s 
radiative-convective boundary which is largely decoupled 
from the nebula. From Figure[4]we infer the following rule of 
thumb: a core must be at least 2 A4® to accrete a volumetri- 
cally significant atmosphere within 1 Myr. (A gas depletion 
timescale shorter than 1 Myr would correspond to a larger 
core mass for the rule of thumb.) For example, a 2 A4® core 
located at 0.3 AU within a disk having p gas = 1600 x 1CP 12 
g/cm 3 accretes a GCR ~ 1% atmosphere within 1 Myr. If 
said core has an Earth-like composition, corresponding to a 
radius of 1.2 R®, a 1% atmosphere would infla te the planet’s 
total radius to 2 R® IjLopez fc Fortnevll2014l l. For context, 
core masses > 2M® also avoid complete atmospheric loss 
by photoevaporation at a ~ 0.3-1 AU, assuming a 10% effi¬ 
ciency in converting incident XUV radiation from the host 
st ar into the kinetic energ y of a pl anetary outflow (F igure 2 
of lLopez fc Fortnevll201 3; see also I Owenfc Wai2013h . 


a (AU) 

T(K) 

-Mcore (-M0) 

Pgas (10 12 g/cm 3 ) 

E 2 ,i (g/cm 2 ) 

0.7 

600 

1 

65 

41 




0.65 




2 

100 

64 




1.0 




4 

160 

102 




1.6 




10 

300 

188 




3 


0.3 

900 

1 

1000 

62 




10 a 




2 

1600 

98 




16 




4 

2600 

156 




26 




10 

4700 

287 




47 



a The {1 AT®, 0.3 AU, low p gas } model could not be evolved because its 
outer radiative zone was too thin to be resolved. The omission of this 
model is unimportant because we do not expect it to yield a volumet- 
rically significant atmosphere, judging from our other 1-M® models. 


We now use the results from the atmospheric accretion sim¬ 
ulations (Figure [4j to identify which of the cores grown in 
the TV-body simulations (Figures [2] and [3]) acquire volumi¬ 
nous gas envelopes vs. stay purely rocky. We carry out two 
methods of analysis. The first and simpler procedure is to 
apply the rule of thumb described in Section ^. 21 We declare 
that all cores that grow to 2 M® within 1 Myr accrete a gas 
envelope, and that all cores that fail to grow that quickly re¬ 
main rocky. We plot the results of this experiment in the top 
panels of Figure [5] for both our Hill-scaled and sound-speed- 
scaled coagulation simulations. We see that only disks with 
higher solid surface densities — E 2 ,i > 50 g/cm 2 — always 
spawn > 2M® cores fast enough to pick up atmospheres. 

The results of a second, more detailed procedure 
are shown in the bottom panels of Figure [5] Here each 
core is prescribed to accrete an atmosphere based on 
the evolutionary tracks shown in Figure [4] By interpo¬ 
lating these tracks, we infer a gas-to-core mass fraction, 
GCR[f, M cor e(f), p S ns(t), a(f)], where GCR[t = 0] = 00 Let 
{tmerger.i} be the array of merger times for a given core. 
For simplicity and because the dependence of the accretion 
rate on the gas density is weak, we assume that the gas disk 
(which has already been depleted to E gas = E 2 based on the 
considerations in Section 2) remains in place for 1 Myr (the 
results are not sensitive to this exact time) and disappears 


2 Technically GCR[t = 0] > 0 because isolation-mass embryos 
can accrete gas before the giant impact era begins. In practice, 
because the gas accretion rate depends so steeply on core mass 
(see Figure 0 , this initial amount of gas acquired is negligible 
compared to that accreted later when the core mass is larger. 


3.3 Combining core growth and atmospheric 
accretion 
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M core =1 M @ M core =2 M @ M core =4 M e M core =10 M 



L Z ,1 (g/cm 2 ) L Zi1 (g/cm 2 ) L z , (g/cm 2 ) L z1 (g/cm 2 ) 

Figure 2. Time to grow to M C o re (as listed on the top of the figure) as a function of solid surface density normalization E Zj i. The results 
shown here are computed from iV-body integrations with embryo masses defined by annuli IORh wide, initialized with random velocities 
uh- The dashed horizontal line marks 1 Myr, an approximate remaining lifetime for the gas disk. The light gray region demarcates those 
surface densities for which M em b = M co re at t = 0 (i.e., no mergers required to grow to M CO re) from a = 1AU (left boundary) to 
a = 0.15 AU (right boundary). In the dark gray region, all embryos at all semi-major axes start with M em b ^ M core and therefore no 
points are plotted. Blue points have 0.5 < a(AU) < 1 and black points have 0.15 < a(AU) < 0.5. In the lower solid surface density disks, 
formation timescales are sometimes much longer than 1 Myr (top left of each panel). 


M core =1 M @ M core =2 M @ M core =4 M @ M core =10 M @ 



Figure 3. Same as Figure |2] but for embryo masses defined in annuli 2c s /Q wide, initialized with random velocities c s . The growth 
timescales exhibit qualitatively the same features as in Figure [2l except they depend less sensitively on semi-major axis (the black and 
blue points nearly overlap); embryo masses defined using c s do not vary as strongly with semi-major axis. 
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t (Myr) 


Figure 4. Gas-to-core mass fractions (GCRs) vs. time for four 
core masses (labeled) at {a = 0.3 AU, T = 900 K} (gray) and 
{a = 0.7 AU, T = 600 K} (black) with high (solid) and low 
(dashed) ambient disk gas density p gas (see Table 1). The red 
striped line illustrates an example “evolutionary track” of a body 
that grows from 1 to 2 to 4A1 q through non-destructive merg¬ 
ers with bodies similar in size to itself; upon each collision, we 
assume that the body loses half its atmosphere and then contin¬ 
ues accreting gas along its new track. By the time the gas disk 
clears at t = 1 Myr (vertical dotted line), our example 4-Mj 
core has accreted a GCR ss 4% atmosphere, enough to make it a 
voluminous mini-Neptune. 


thereafter: 


Pgas — E z ll/Cs 
= 0 


for t < lMyr 
for t > lMyr . 


If the first merger f m erger,i > 1 Myr, the final atmospheric 
gas content of the planet equals GCR[t = 1 Myr, M core = 
M e mb]- Otherwise we interpolate the evolutionary tracks to 
calculate GCR[f = lmer g er,i] just before the first merger. 
Just after the first merger, we assume that, because of 
atmospheric loss from the giant impact, the new bod 
starts with half the GCR of the larger merging body[ 
This is meant to be an app roximation to the more realis¬ 
tic tre atments developed b y ISchlichting. Sari fc Yalinewicbl 
J2015I) and llnamdar fe Schlichtind ( 20151 ). Note that the 
overwhelming majority of mergers are between bodies of 
comparable mass. From there, the GCR follows a new track 


3 In practice, this atmospheric loss is not important because gas 
accretion occurring after the last merger prior to the disappear¬ 
ance of the gas disk typically dominates. 


appropriate to its new M core , accreting additional gas for 
t = min(t me rger, 2 ,1 Myr) — f merger, 1 - This process repeats 
until fmerger,* > 1 Myr. See Figure 4 for an example evolu¬ 
tionary trajectory. 

Planets that attain GCR is 1% are considered gas- 
enveloped and colored orange in Figures [5] and [6] This sec¬ 
ond method yields essentially identical results to the first 
rule-of-thumb treatment. 

In all panels of Figure[ol we see that the lowest E 2 disks 
produce predominantly rocky (black) planets; intermediate 
E z disks yield a combination of gas-free and gas-enveloped 
(orange) planets; and the highest E 2 disks spawn exclusively 
gas-enveloped planets. Figure [5] indicates that virtually all 
planets with final masses above ~5M® accrete gaseous en¬ 
velopes, suggesting that radial-velocity discoveries like GJ 
87 6 e (a 12. 5 Mm planet with a 124-day orbital period; 

I Rivera et al.l I 2 OI 0 I ) or 61 Vir c (a 10.6 M® planet with a 
38-day orbital period; lYoet et al.l l2Q10h are mini-Neptunes 
rather_than super-Earths. 

llnamdar fe Schlichtind (120151 ) also recently simulated 
the atmospheric mass fractions attained by an ensemble of 
planets that grow from isolation mass to core mass during 
the giant impact stage. Their isolation masses are formed 
from annuli 5J?h wide and having a solid surface density E 2 
such that observed exoplanets of mass M v can form from 
feeding zones tha t are A a = 2v esc /Q, wide, where u esc is 
evaluated for M p (ISchlichtindl2014l ). They simulate the ac¬ 
cretion of gas during the isolation mass stage; subsequent 
atmospheric loss during giant impacts; and the final post¬ 
giant impact accretion of gas. The giant impacts are based 
on Monte Carlo simulations — by comparison with the N- 
body simulations presented here — and atmospheric loss is 
treated in a more detailed and realistic way than our sim¬ 
ple factor-of-2 prescription. When they consider only atmo¬ 
spheres accreted during the isolation-mass stage and subse¬ 
quently eroded through impacts, the resulting atmospheric 
masses fractions are low, on the order of ~10 -3 (their Fig¬ 
ure 8). However, atmospheric mass fractions reach several 
percent for subsequent accretion by core masses above 2 
beyond 0.15 AU (their Figure 9), unless the energy from 
giant impacts happens to be released over exactly the disk 
dissipation timescale. Thus our results are consistent with 
their post-giant impact results, in which gas accretion is 
dominated by the largest core mass achieved before the dis¬ 
sipation of the gas disk. 


3.4 Effects of migration 

The formation theory we have explored above is an in-situ 
theory; the planets are assumed to form in place. We sketch 
here how our results are impacted by orbital migration. 
Transporting bo dies by gravitat ional torques in a disk (for 
a review, see Kiev fe Nelsonll2012l ) can: (1) supply planetary 
embryos from the outer disk to the inner disk; (2) estab- 
lish a resonant c hain of protop lanets in the inner disk (but 
see iGoldreich fe Schlichtinell2014l for reasons to believe such 
resonant locks are easily broken); and (3) move fully formed 
super-Earths and mini-Neptunes to the inner disk. We will 
argue that all these effects are compatible with our thesis 
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♦ rocky » gas-enveloped 



Figure 5. Compositional outcomes and how they depend on solid 
surface density (as measured by its normalization S 2j i evaluated 
at 1 AU). Top panels are based on our “rule-of-thumb” treatment 
for gas accretion: a planet is considered gas-enveloped (colored or¬ 
ange) if it grows to 2in < 1 Myr; otherwise it stays purely 
rocky (black). Bottom panels are derived from the detailed evo¬ 
lutionary tracks shown in Figure |4] planets are considered gas- 
enveloped if they attain a GCR ^ 1% within 1 Myr; see text 
for details. Left panels utilize data from the TV-body simulations 
shown in Figure [2] which are based on Hill scalings; right panels 
derive from Figure |3] which uses sound-speed scalings. How the 
compositional outcome depends on solid surface density is insen¬ 
sitive to modeling specifics, as judged by the similarity of all four 
panels. 

that the disk’s solid surface density is a major determinant 
in a planet’s final gas-to-rock ratio. 

In scenario (1), planet formation would proceed as de¬ 
tailed above (Sections 3.1-3.3) but with E 2 reflecting the 
surface density of embryos supplied from the outer disk 
rather than the surface density of dust/planetesimals in the 
inner disk. The simulations of Sections 3.1-3.3 model events 
post-dating the growth or delivery of embryos, so migra¬ 
tion scenario (1) would not impact our conclusions quali¬ 
tatively. Scenario (2) would be characterized by larger or¬ 
bital spacings between protocores and longer orbit crossing 
times. Nevertheless the de-stabilization time should depend 


on planet masses (see, e.g., t he resonant chains exhibited 
by th e inner moons of Uranus; iFrench, Dawson Sz Showalterl 
l2015h which in turn would reflect the isolation masses in 
the outer disk; thus merger times should still depend sen¬ 
sitively on E 2 ,i. Finally, under scenario (3), coagulation of 
rocky cores and accretion of gas would occur at larger or¬ 
bital distances, where isolation masses are larger and or¬ 
bital timescales are longer. The latter two effects appear to 
cancel for our fiducial semi-major axis range (0.04-1 AU); 
the black vs. blue points in Figures 2 and 3 exhibit similar 
trends. Whether they cancel so nearly at larger semi-major 
axes remains to be determined. If they do not cancel, then 
the threshold E 2j i required to produce gas-enveloped plan¬ 
ets could be larger or smaller, but it would still exist; i.e., 
we would still expect a population of purely rocky planets 
in the lowest E 2 p disks. Because feeding zone annuli tend 
to be wider at larger orbital distances, cores formed in the 
outer disk do not need as large a E z> i to produce cores that 
are massive enough to accrete significant amounts of gas as 
isolation masses. 

Recent models incorporating the effects of migration 
have assessed the effect of a parameter analogous to 
S z .i on t he pr opert ies of super-Earths and mini-Neptunes. 
ICossou et al.l (2014) simulated planet formation with N- 
body integrations combined with prescriptive migration 
maps. Not treating gas accretion by cores, they found that 
disks with a higher mass in embryos delivered more mas¬ 
sive super-Earth cores to the inner disk. We would expect 
these more massive cores to more eas ily acquire atmospheres 
(Section 3.2). In addition. Ilda fe Linl (l2010l l found that, for a 
range of migration speeds, the typical mass of super-Earths 
that migrate to < 1 AU increases with disk solid surface den¬ 
sity normalization. The delivery of higher-mass planets to 
the interior regions of disks with overall higher solid density 
is consistent with our statement that the masses of planets 
that arrived by migration reflect the isolation masses in the 
outer disk. 


4 COMPARISON TO THE KEPLER SAMPLE 

We argued above that the disk’s solid surface density is crit¬ 
ical for determining whether a planet forms rocky or gas- 
enveloped. Here we translate the results from Section 3 — 
compositions as a function of E 2 - - into observables to com¬ 
pare with the large Kepler sample of super-Earths and mini- 
Neptunes. 

4.1 Observational proxies 

Kepler has discovered about 3000 planet candidates with 
sizes between that of Earth and Neptune. Statistical mod¬ 
eling of the potential astrophysical false positive population 
has reveal ed that less than 10% of the candidates are false 
positives llMorton fe Johnsonll20lil : iFressin et al .1 [20131 ). al¬ 
lowing us to treat the sample as representative of true plan¬ 
ets. Most Kep ler candidates have me asured radii but un¬ 
known masses. [Lopez fe Fortnevl(l2014l l recently argued that 
a planet’s radius can serve as a proxy for its composition 
(i.e., gas-to-rock content), assuming an envelope composed 
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of hydrogen, helium, and an admixture of metals@ In par¬ 
ticular, planets larger than 27?® require a gas-to-core mass 
fraction GCR exceeding ~1%. This prediction has been 
borne out by the Kepler data. From a stat istical analysis 
of plane ts with measured ma sses and radii, Rogerj (2015) 
(see also I Weiss &; Marcvll2014l and lWolfgang &; Loped 2015 1 
found that the cut-off radius between planets dense enough 
to be rocky and those requiring a significant volatile compo¬ 
nent sits at ^1.67?®; a planet of radius 1.67?® has a ~50% 
probability of being purely rocky. Between 1.5 to 2 7?®, plan¬ 
ets appear to transition from being purely rocky to having 
volum etrically significant gas envelopes (Figure 5 of [Roger j 
l2015ll . A similar transition is suggested by planets with or¬ 
bital periods shorter than 1 day — a.k.a. ultra-short period 
planets or USPs. Because USPs are so intensely irradiated 
by their host stars, they may have lost their primordial at¬ 
mospheres to photoevaporation and may therefore exhibit 
a radius distribution correspon ding to that of purely rocky 
cores. All the USPs studied by ISanchis-Oieda et all (|2014T ) 
have radii 7? < 1.687?®, consistent with the interpretation 
that planets with larger radii possess voluminous gaseous 
envelopes. 

Together, these recent advances allow us to identify Ke¬ 
pler planets as rocky or gas-enveloped based on their radii 
alone. We pay particular attention below to Kepler plan¬ 
ets with periods longer than ~15 days, since photoevapo- 
ratio n can remove the envel opes of shorter period planets 
(e.g., iLopez fe Fortnevll2013l ). leaving behind a planet that 
is purely rocky by nurture rather than by nature. 

Having opted for planet radius as a proxy for 
composition, we now need a proxy for E z . Motivated 
by the strong correlation bet ween giant pla net occur¬ 
rence and host star metallicit y (lSiMitos J Jhraelian_fchlayoj 
200 ll: Fischer fe Valentil [20051: IBuchhave et al.l 120121 . 1201? 


Wang fe Fischeil 20151 ). we use the spectroscopic [M/H] 


of planet-hosting stars as a proxy for the surface density 
E z of their primordial planet-forming disks. Kepler follow¬ 
up teams have observed stellar spectra for h undreds of 
KOI (Kepler Object of Interest) host stars. IBuchhave et al.l 
(I2014T) recently published a catalogue of [M/H] measured 
using t he stellar parameter cla ssification (SPC) tool devel¬ 
oped by IBuchhave et all d2012l) . Because the host stars span 
a narrow stellar mass range, while [M/H] spans nearly a full 
dex, we expect that [M/H] should track E z — or at least 
its value averaged over the entire disk (it may not track the 
local value of E z because solids can concentrate radially in 
disks, so that the local disk metallicity can deviate from the 
bulk stellar met a llicity — more on this point in Section 4.3). 
IBuchhave et al.l (l201df ') reported that hosts of KOIs with 2- 
4 7?® have higher metallicities than hosts of smaller KOIs. 
ISchlaufmanl J2015T ) recently argued against the statistical 
significance of this finding. 


4 In many models, the atmospheric metallicity can be as large as 
Z ss 0.4. 


♦ rocky « gas-enveloped 



Figure 6. Compositions of planets (orange = gas-enveloped, 
black = purely rocky) spawned by different solid surface den¬ 
sity disks across a range of orbital periods. Top: based on the 
results shown in the bottom left panel (Hill scalings + detailed 
gas accretion prescription) of Figure[5] Bottom: based on the re¬ 
sults shown in the bottom right panel (sound-speed scalings + 
detailed gas accretion prescription) of Figure [5] 


4.2 A metallicity divide for rocky planets 

In Figure [6] we plot the results of the simulations from Sec¬ 
tion 3 with the added dimension of orbital period. The tran¬ 
sition E z that divides gas-enveloped planets (orange points) 
from purely rocky planets (black points) decreases slightly 
with orbital period. Longer orbital periods correspond to 
larger isolation masses, lowering the threshold E z to spawn 
gas-enveloped planets only. This effect may be related to 
the observed rise in the occurr ence rate of Neptu ne-sized 
planets at longer orbital periods dDong fc Zhul2013l ). though 
the concomitant weakening of photoevaporation at large or- 
bi tal distances l i kely also cont ribut es. (S ee also Figure 7 
of lForeman-Mackev, Hogg fe Mortonl l2014l.) Regardless, be¬ 
yond 20 days, the transition E Z: i flattens to about 40 g/cm 2 , 
about 2x the lowest E z ,i that can spawn close-in planets 
more massive than Earth, or about 4x the value of the 
minimum-mass solar nebula. Note that near this transition 
E Zj i, a mixture of gas-enveloped and purely rocky planets is 
produced. 

Next we plot, in Figure 0 the observed host star 
metal fraction (r elativ e to solar) vs. orbital period using the 
IBuchhave et al.l (120141 ) catalogue of stars with spectroscopic 
parameters. We use radius as a proxy for composition, color¬ 
coding orange those planets with radii exceeding 2 7?® to 
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P (days) 


Figure 7. Metal-rich stars lack purely rocky (R < 1.577® = 
black) planets on long er period orbits (P > 15 days); data from 
IBuchhave et ali 1120141 ). We focus on orbital periods > 15 days 
to mitigate the influence of photoevaporation. Metal-rich stars 
with predominantly gas-enveloped planets (pink box) can be dis¬ 
tinguished from metal-poor stars with a combination of gas- 
enveloped and purely rocky planets (gray box). We argue that 
higher metallicity stars are accompanied by disks with higher 
solid surface densities, which in turn spawn > 2Mj cores faster, 
within the gas disk lifetime of ~1 Myr; these planets more readily 
acquire atmospheres and inflate their radii to > 277® (orange). 
The dotted line indicates solar metallicity for reference. 


signify that they likely have volumetrically significant at¬ 
mospheres. We remove KOIs designated as false positives 
in the NExSci database dMullallv et al.ll2015l) . Figure [7| at¬ 
tests that at orbital periods > 15 days, beyond the reach 
of photoevaporation, metal-rich stars (in the pink box) lack 
rocky planets, whereas lower metallicity stars (in the gray 
box) host a mixture of rocky and gas-enveloped pla nets. 
This trend is related to that reported by IBuchhave et al.l 
ll20l4) . who found that planets above ~1.7 Earth radii orbit 
higher metallicity stars. Instead of identifying a cut in planet 
radius, we identify a cut in metallicity above which stars 
host exclusively gas-enveloped, not rocky, planets. Unlike 
the simulations shown in Figure 6, the observations shown 
in Figure 7 show no evidence that only rocky planets orbit 
the lowest metallicity stars; we will expand upon this point 
in Section 4.3 (see Figure 8). 

We compute a K-S statistic comparing the distributions 
of planetary radii in the pink vs. gray boxes in Figure [7] 
yielding a formal statistical probability of 0.00045 that the 
difference in the distribution of radii is due to chance. It is 
not clear theoretically what the exact value of the metallic¬ 
ity divide Z^iv should be between the pink and gray boxes; 
the existence of a metallicity divide is physically motivated 
by Sections 2-3 which highlight the role played by the solid 


surface density £ z in determining a planet’s final gas con¬ 
tent, but the exact mapping between £ z and [M/H] is not 
known. Despite this uncertainty, the high (formal) statistical 
confidence reflected in the K-S test is robust against a range 
of values for Zdiv- In the Appendix, we perform an alterna¬ 
tive statistical test that accounts for our freedom in choosing 
Zdw; we find there that a model with a metallicity divide is 
preferred to one without at 95% confidence. The period cut 
of 15 days delineating the left boundaries of the pink and 
gray boxes is motivated by considerations of photoevapo¬ 
ration; for example, photoevaporation has been invoked to 
explain the rocky comp osition of Kepler-36b, which has a 
13.8 day orbital period llLopez fc Fortnevll2013l ). 

In interpreting the observations, we have assumed that 
the reported metallicities are correct relative to one another. 
Figure 7 should be remade if improved stellar parameters are 
obtained. For example, if the stars reported as metal-rich 
were actually giant stars, it would not have been possible 
to detect small planets orbiting them, and their identifica¬ 
tion as planet-hosting stars would be incorrect. Another po¬ 
tential complication is that stars more rich in metals tend 
to be larger than their metal-poor counterparts, making it 
more difficult to detect smaller (rocky) planets orbiting the 
former. In the sample here, the average size of stars that 
have metallicities Z* > 1.5Z® and that host a 1-4 7?® KOI 
with a period > 15 days is 1.3 ± 0.477®. By comparison, 
for stars with Z < 1.5Z®, the average size is 1.1 ± 0.477®. 
Nevertheless, the stellar radius ranges of the two samples 
have substantial overlap, and the difference in the median 
stellar projected areas implies only a ~20% difference in 
signal-to-noise for transit detections. We also tried restrict¬ 
ing the sample to stars for which otherwise identical planets 
of 1.377® could have been detected; this reduced the sample 
size but did not change the results qualitatively. Therefore 
we interpret the trend identified in Figure [3 to be real and 
not just a selection effect, but we recommend revisiting the 
trend as the sample size of spectroscopically measured KOI 
hosts grows. We also recommend the publication of the joint 
two-dimensional posteriors for stellar metallicity and stellar 
radius so that any covariance in the inferred parameters can 
be accommodated. 


4.3 Planet radius vs. host star metallicity in detail 

We can add another layer of modeling by computing precise 
planetary radii from our simulation s. The interior structure 
models of iLopez fe Fortnevl (|2013l ) yield an empirical for¬ 
mula for planetary radius as a function of core mass and 
GCR. Planet radii so evaluated using the simulation data 
from Section 3 are plotted against E Zj i in Figure [S] and com¬ 
pared alongside the observed distribution of planet radii vs. 
host star metallicity. In both simulations and observations 
we see an “empty wedge” at large £ z (large host star metal¬ 
licity) and small planetary radius. 

However, Figure [8] also reveals a discrepancy with 
the observations. The simulations predict a second empty 
wedge: an absence of large radius planets in the lowest £ z 
disks. But this region of parameter space is actually pop- 
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Figure 8 . Planetary radius as a function of solid surface density (its normalization S 2j i at 1 AU) and its observational proxy, the 
host star metallicity (Z/Zq normalized to solar). Middle and right panels are take n from t he simulatio ns underlying the bottom panels 
of Figure 5 (left and right panels, respectively), combined with the models of lLopez fc Fortney! ||2013| V Blue wedges are mostly empty 
regions of parameter space. Theory and observation agree on the lower right wedge — the absence of small radius (i.e, purely rocky) 
planets orbiting stars of high metallicity (i.e., born in disks of high solid surface density). But the upper left wedge nominally predicted 
by theory is actually filled by observations — perhaps because radial drift of solids in the primordial disk renders stellar metallicity an 
inadequate proxy of disk solid surface density, or because the upper left wedge is populated by planets that migrated from larger orbital 
distances. 


ulated by the observation^ There are a few possible ex¬ 
planations. First, [M/H] may be an imperfect proxy for 
E,,i. Some disks orbiting metal-poor stars may have dis¬ 
proportionately large E 2 if the solids in their primordial 
disks drifted radially inward by aerodynamic drag and ac¬ 
cumu l ated at small orbi t al dist a nces (e.g., Youdrn&Shu| 
200^_[Youdin_fciCdnjmgl__|200^JChja^ 120101 : 


Hansen & Murray boi^ r ichatteriee fe Tanl 12014 1. In other 
words, the disk’s solids-to-gas ratio at any particular loca¬ 
tion can differ from the stellar metalli city because solids an d 
gas can segregate in disks (see also [Andrews et al.l 120121) . 
Second, gas-enveloped planets may have formed at larger 
separations beyond 1 AU — where isolation masses are 
larger — and migrated in (see Section 3.4). 


5 CONCLUSIONS 

Kepler has discovered an abundance of planets with volu¬ 
minous atmospheres. A small mass per centage of volatiles 
— 1% for typical mini-Neptunes (e.g., IWolfgang fe Lopea 
120151) — can dramatically inflate a planet’s radius beyond 
that of its rocky core, with significant consequences for sur¬ 
face temperature, pressure, and habitability. Gas-enveloped 
planets are the rule around stars with supersolar metallici- 
ties; atmosphere-laden planets are also found orbiting stars 
with subsolar metallicities, together with purely rocky plan¬ 
ets having practically no contribution to their radii from 
their atmospheres. 

What determines whether a rocky core acquires an at¬ 
mosphere or not? We have found that the protoplanetary 
disk’s surface density in solids is a key ingredient in recipes 


5 ISchlaufmanl (' 201 oh found a linear correlation between planet 
radius and host star metallicity using a collection of candidates 
with sizes extending to 15 ft®; that trend may be dominated by 
planets larger than 4 ft® so we do not compare here. 


for forming rocky vs. gas-enveloped planets. In Section 2 we 
used order-of-magnitude scaling relations to show that low 
solid surface densities prolong the timescale for a core to 
grow from embryos, so much so that by the time the core is 
massive enough to acquire an atmosphere, there may not be 
any disk gas left. Long coagulation timescales in low solid 
surface density disks are a consequence of the larger number 
of mergers required to assemble a core from smaller isolation- 
mass embryos. In Section 3, we combined A-body simula¬ 
tions of core assembly with one-dimensional gas accretion 
models to show that the highest surface density disks pro¬ 
duce primarily gas-enveloped planets, whereas lower solid 
surface density disks can produce b oth gas-enveloped and 
rocky planets. In Section 4, using the lBuchhave et al.l (|2014l ) 
spectroscopic sample, we presented observational evidence 
for a lack of purely rocky planets orbiting metal-rich stars 
at orbital periods > 15 days (Figure 7 and left panel of 
Figure 8). This trend is another manifestation of the cor¬ 
relation between planet radius and host star metallicity 
for small planets repo rted previously jBuchhave et al.|l2014l ; 

I Wang &T~F ischcr 2015). If we assume a 1-to-l correlation be¬ 
tween host star metallicity and disk solid surface density, 
then we can reproduce the absence of purely rocky plan¬ 
ets around metal-rich stars (i.e., high solid surface density 
disks; Figure 8). But this same simple model does not re¬ 
produce the observation that metal-poor stars also host gas- 
enveloped planets with sizes up to 4 Earth radii (Figure 
8). This probably means that the assumption of a 1-to-l 
correlation between stellar metallicity and disk surface den¬ 
sity is inadequate; radial drift and accum ulation of solids 
(for a review, see IChiang fc Youdinl l201tf > can yield high 
surface density disks (and therefore gas-enveloped planets) 
even around stars with low bulk metallicity. Alternatively, it 
might also be that the in-situ planet formation models used 
here, which have a hard time producing “super-puffy” plan- 
ets wit h low core masses and extended atmospheres (e.g., 
[m asu as, need to account for migration. As argued in 
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Section 3.4, migration should produce qualitatively similar 
trends in planet composition vs. disk solid surface density. 
The subset of mini-Neptunes orbiting metal-poor stars may 
have migrated from afar, where isolation masses are larger. 

These results may bear on other recently discovered 
trends between host star metallic i ty an d the properties 
of small planets. lAdibekvan et all l|2013l . their Figure 1) 
found that low-mass planets (< 10M®) orbiting metal- 
rich stars have orbital periods shorter than ~20 days, 
whereas more massive planets, and planets orbiting metal- 
poor stars, span a range of orbital periods. The restric¬ 
tion of mini-Neptunes/super-Earths to short orbital periods 
around metal-rich stars may be a manifestation of how em¬ 
bryo masses — and presumably final core masses — increase 
with increasing orbital distance. That is, at large orbital dis¬ 
tances in a high surface density disk, embryo masses may 
be so large that the disk canno t help but spawn hig h -mass 
(> 10M®) planets. Relatedly, iBeauge fe Nesvornfl (|2013l . 
their Figure 6) found that super-Earths and mini-Neptunes 
(radii < 47?®) hosted by low-metallicity stars have orbital 
periods greater than 5 days; that super-Earths and mini- 
Neptunes with orbital periods shorter than 5 days are hosted 
by high-metallicity stars; and that metal-poor stars do not 
harbor planets with radii > 47?®. Their reported trends 
— which may stem from the higher solid surface densities 
expected to accompany higher metallicity stars, in concert 
with the rise in embryo mass with orbital period — could be 
revis ited with the latest sam ple of spectroscopic metallicities 
from iBuchhave et al.l (l2014fl . 

We have assumed in this paper that mini-Neptunes are 
either purely rocky or gas-enveloped. The possibility that 
mini-Neptunes are “water-worlds” having a large fraction of 
their mass in water or ices cannot be discounted based on 
mass and radius measurements a lone (see, e.g., the case of 
GJ 1214b: [Rogers &; Seageill2010lh Our study favors the in¬ 
terpretation that mini-Neptunes are gas-enveloped and not 
water-worlds; we would not expect the prevalence of purely 
rocky vs. icy planets to be a strong function of disk solid 
surface density or host star metallicity. 

Our study supports the idea that the volatiles in mini- 
Neptunes are directly accreted from the primordial nebula 
rather than outgassed from rock, as outgassing does not 
obviously lead to the observed trends with stellar metal¬ 
licity / disk solid surface density. Further evidence for 
nebular accretion can be fou nd in the ultra-short perio d 
(< 1 day) planets studied by ISanchis-Oieda et al.l (1201414 . 
These have a radius cut-off of ~1.7 7?® — consistent 
with th e tran sition from rocky to gas-enveloped planets 
dRogersI [2015l l. but below the ~3 7?® break above which 
the pla net occurrence rate decreases at larger orbital dis¬ 
tances (IPetigura. Marcv fc Howardl l2013h . The absence of 
planets with voluminous gas envelopes at the shortest or¬ 
bital periods is more consistent with nebular accretion: at 
these close-in distances, primordially accreted atmospheres 
are readily lost to photoevaporation d uring the firs t ~0.1 
Gyr of the host star’s evolution (e.g., lOwen fc Wull2013lh 
whereas steam or gradually outgassed atmospheres can be 
maintained or replenished throughout the star’s lifetime. 
The one notable exception of a non-rocky ultra-short period 
planet is 55 Cnc e, which has an orbital period of 0.74 days 


dPawson fc Fabrvckvll2ofoi ), a radius of 2 7?® , and a density 
too low to be purely rocky (IWinn et al.M 201 14. Intriguingly, 
it is situated in a system with two close-in giant planets 
whos e proximity to a mean-motion resonance dMarcv et al.l 
120021 : iMcArthur et al.l l2004 iFischer et al.l l2008h suggests 
they underwent orbital migration. Thus 55 Cnc e may also 
have been transported from afar; it may be one of an un¬ 
common class of migrated icy planets, interloping among 
the majority of super-Earths and mini-Neptunes that coag¬ 
ulated and acquired their atmospheres in situ. 
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APPENDIX A: STATISTICAL SIGNIFICANCE 
OF ABSENCE OF ROCKY PLANETS 
ORBITING METAL-RICH STARS 

As an alternative statistical test that accounts for our free¬ 
dom in choosing the metallicity divide, we compute an odds 
ratio to compare a model with a metallicity divide (model 
1) to one without (model 2). For the former, we marginalize 
over the metallicity divide. The odds ratio of model 1 over 
model 2 is: 

Pmodell|data P<lata|model 1 ^ Pmodell (Al) 

Pmodel2|data /f J ala | model 2 Pmodel2 

If we give equal prior weight to the two models, the odds 
ratio becomes: 

Pmodelljdata Pdata| model I , , .. , 

- - . 

Ptrioilel 21 data Pdata|model2 

Here the data are the inferred compositions of observed 
planets, so we define a probability that an observed planet is 
rocky, p r ock y , given its observed radius R p ,i and uncertainty 
o\R P ,i. The uncertainty OR V ,i is assumed to define a normal 
distribution and is comp uted by propagating th e uncertainty 
from the stellar radius (iBuchhave et al.l 12014! ) and transit 
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depth dMullallv et al.ll2015l ~). assuming that the uncertain¬ 
ties in these quantities are independent. For simplicity, we 
estimate 


Prockylflp — 0, Rp > 27?® 
1, R p < 277® . 


(A3) 


W e clarify t hat this p ro ck y differs from the p roc k y defined 
by iRogera d2015l l. who compute p r0 cky based on both mass 
and ra dius measurem ents. Our p roc ky abstracts the re¬ 
sults o fiRogersI d2015l'l and other st udies (| Lopez fc Fortnevl 

"" 1 


2014 


2014 S anchis-Oieda et all l20l4 IWeiss fc Marcvl 
Wolfg ang &: Lopez 20151 1 to assert that 27?® is an approx¬ 


imate cut-off above which a planet is very unlikely to be 
purely rocky. We marginalize over the uncertainty in the 
planet’s radius: 


Procky 




exp 


(7?. — R P ,i ) 2 \ 


dR. (A4) 


Next we evaluate the denominator of the odds ratio in 
Equation IA21 the probability of the data given the model 
without the metallicity divide (model 2): 


Pdata|model2 

lit / [/rockyProcky,* ~\~ (1 /rocky) (1 Procky,*)] dfrocky (A.5) 

Jo 

where i refers to an individual planet. Equation IA5I 
marginalizes over /rocky, the fraction of the sample com¬ 
prised of rocky planets (e.g., / roc ky = 0.5 for an even mixture 
of rocky and gaseous planets). 

A similar formula to Eauation lA5l vields the probability 
of the data given model 1 for host stars with metallicities 
below the divide Zdiv', now / r0 ck y is the fraction of the sample 
below the divide that is rocky: 

Pdata^ |modell,2* > i<Z’ c ji v 

I [/rockyProcky ,* T (1 /rocky) (1 Procky,*)] d/ roc ky . (A6) 

JO 

By contrast, above Zdiv, planets are unlikely to be 
rocky. We marginalize over /r OC ky! the fraction of the sample 
above the cut that is rocky, assuming f ' acky ranges from 0 
to 10%: 


Pdata^ |modell,.Z* ) i>^di v 

1 

[/rockyProcky,* (1 /rocky) (1 Procky,*)] c/rocky • (A7) 

Marginalizing over Zdi v , we have 



iii (p, 


'dataj |modell,Z* 


Pdata| modell 



iv Pdataj |modell,Z* ; j >Zdiv ) • (A8) 


We use the observed sample (Figure [7]) to evaluate 
Equations IA5I and IA8I We thereby obtain an odds ratio 
(Equation IA2I) of 20 in favor of the two-population model 
(i.e., 95% confidence). The integrand in Equation IA8I peaks 
for Zdw = 1.5. 
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